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Abstract. An analysis of ages, metallicities, and [a/Fe] ratios of globular cluster systems in early-type galaxies 
is presented, based on Lick index measurements summarized in Puzia et al. (2004, Paper I of this series). In the 
light of calibration and measurement uncertainties, age-metallicity degeneracy, and the relative dynamic range of 
Lick indices, as well as systematics introduced by abundance ratio variations (in particular variations of [a/Fe] 
ratios), we find that the most reliable age indicator for our dataset is a combination of the Lick Balmer-line 
indices H7A, H/3, and US a- [MgFe]' is used as a spectroscopic metallicity indicator which is least affected by 
[a/Fe] variations. We introduce an interpolation routine to simultaneously derive ages, metallicities, and [a/Fe] 
ratios from diagnostic grids constructed from Lick indices. From a comparison of high-quality data with SSP 
model predictions, we find that ~2/3 of the globular clusters in early-type galaxies are older than 10 Gyr, up to 
1/3 have ages in the range ~ 5—10 Gyr, and only a few cluster are younger than ~5 Gyr. Our sample of globular 
clusters covers metallicities from [Z/H] w —1.3 up to ~ 0.5 dex. We find that metal-rich globular clusters show 
on average a smaller mean age and a larger age scatter than their metal-poor counterparts. [a/Fe] diagnostic 
plots show that globular cluster systems in early- type galaxies have super-solar a/Fe abundance ratios with a 
mean [a/Fe] = 0.47 ± 0.06 dex and a dispersion of ~ 0.3 dex. We find evidence for a correlation between [a/Fe] 
and metallicity, in the sense that more metal-rich clusters exhibit lower a-element enhancements. A discussion 
of systematics related to the Lick index system shows that the method suffers to some extent from uncertainties 
due to unknown horizontal branch morphologies at high metallicities. However, these systematics still allow us 
to make good qualitative statements. A detailed investigation of indices as a function of data quality reveals 
that the scatter in Balmer index values decreases for higher-quality data. In particular, extremely low Balmer 
index values that are lower than any SSP model prediction tend to disappear. Furthermore, we find that observed 
photometric colors are in good agreement with computed SSP colors using ages and metallicities as derived from 
the spectroscopic line indices. 
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1. Introduction 

Two prominent models of early-type galaxy formation 
are lively discusse d in the literature. In the monolithic- 
collavse scenario llTinslevlll972t lLarsorJ Il975t Isilk1ll977t 
lArimoto fc Yoshiil 19871 etc.) the majority of stars in 
early- type galaxies forms early, at redshifts z <; 2. 
Empirical scaling laws, such as the fundamental plane, 
the Mg-cr relation, and the [ev/Fe]-<7 relation support 
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such an early formation epoch llBower et alJI 


19921 iBcndcr 


1996 


: iBender et alJ 


1996: iTreu et all 1999: 


iKuntschnei 


200C 
2003 


: iTraeer et all 2000albl iThomas et alJl2005l IPeebles 


) . On the other hand, the hierarchical merging picture 


(White & Recs 1978; Kauffmann et al. 1993; Baugh et al. 



19981 ICole et al.ll200f3 ISomerville etal1l200l[) sees early 



type galaxies as the result of multiple merging and ac- 
cretion events of smaller units over an extended pe- 
riod of time until the very recent past. In this way 
a significant fraction of stars is formed below a red- 



tion events (e.g. 


Whitmore & Schwcizerlll995Hlbata 


et al. 


1199,4 iKoo et al. 


19961) and the existence of kinematically 


decoupled cores 


(Franx & Illmgworthll 19881 Bender 




Jedrzeicwski & Schechterl Il988t ISurma & Benderl 




Davies et al.l2001) arc clear lv enforcing' arguments for this 



scenario. 

These models describe two antipodal paradigms of 
galaxy formation. The quantification of their importance 
as a function of rcdshift, environment, and galaxy mor- 
phology is necessary to provide a detailed insight in galaxy 
formation and assembly. A major difference between the 
two pictures are the different star-formation histories of 
early-type galaxies. While in the hierarchical picture most 
massive galaxies are thought to experience long assembly 
time scales, the monolithic collapse scenario predicts very 
early and short bursts of star formation. Clearly, the pre- 
dicted star-formation histories stand in marked contrast 
and are an important piece of evidence to differentiate 
between these two models. 

However, to recover star-formation histories, one ide- 
ally has to resolve the underlying stellar populations which 
build up a galaxy. Only a few early-type systems are close 
enough so th at their halos can be r esolved int o single stars 
(e.g. in M3 2:lGrillmair et al. Il99d NGC 3115 :lElsonll997t 
NGC 5128: lHarris fc Harrisl20oT NGC 3379: lGregg et alJ 
Most photometric and spectroscopic studies of the 
unresolved diffuse light are hampered by the mix of ages, 
metallicities, and abundance ratios in the stellar popu- 
lations (e.g. iMaraston fc Thomas! E?) 00L In combination 
with the well-known age-metallicity degeneracy I Faberl 
Il972t IO'ConneIilll976t IWorthevlll994h . it is extremely dif- 
ficult to disentangle even the different major stellar pop- 
ulations from studies of the diffuse light only, let alone to 
reconstruct a detailed star formation history. 



1.1. Globular Cluster Systems as Diagnostics of 
Star-Formation Histories 

Globular cluster systems are very useful tools to study star 
formation histories of galaxies. Several arguments support 
the hypothesis that globular clusters trace major star- 
formation events in galaxies. (1) The formation of massive 
star clusters, which are likely to survive a Hubble time 
as globular clusters, is obs erved in interacting /merging 
and starburst ga l axies (e.g. IWhitmore fc Schweizerlll99,4 
ISchweizerl Il997t I.Tohnson et al] Il99flt lOoudfrooii et al] 



l200ll life meier et al.ll2002|) . (2) Young massive star clus- 
ters are observed in "simmering" late-type galaxies. Their 
number is correlated with the star formation rate per 
unit area in these systems l(Larsen fc RichtAerl E'oOO'). (3) 
Normalizing the number of globular clusters to the total 
baryo nic mass of the hos t, reveals a surprisingly constant 
value 4McLaughlini ri999) and points towards a tight link 
between star and globular cluster formation. 

In other words, globular clusters are fossil records 
of major star-formation episodes in galaxies. Their ages, 
metallicities, and abundance ratios can provide detailed 
information on the formation histories of their host sys- 
tems. As simple stellar populations, globular clusters con- 
sist of stars characterized by one age and one metallicity. 
Hence, the interpretation of their observed colours and 
spectra is less ambiguous than for the diffuse light. Their 
ubiquity in all galaxy types and their high surface bright- 
ness make them easy to observe out to large distances. 

1.2. Spectroscopy of Globular Cluster Systems 

Spectroscopy opens an independent way to access ages 
and metallicities of globular clusters besides photometry. 
Low resolution spectroscopy (R 1000) provides infor- 
mation on the strength of prominent diagnostic features, 
such as the Balmer series and some relatively strong iron 
and other metal features. Moreover, it can provide clues 
on the basic chemistry of globular cluster systems. 

The Lick index system standardizes the measurement 
of spectrosco pic line indices for many strong absorp- 
tion features llBnrstein et^Llll984 IWorthev et al .111 994 
IWorthev fc Ottavianilll997tlTrager et al.lll998h . The out- 
standing role of the Lick system is its provision of in- 
dex measurements for many stars with a wide range in 
log<7,T e ff, and [Fe/H]. With this information, so-called fit- 
ting functions 1 are computed which are an essential ingre- 
dient of theoretical S SP model predictions . Using index 
response functions toipicco fc Belli Il995l) recent SSP 
models l|Trager et alJ l2000al IThomas et al 1 120031 12004^ 
provide also information on fundamental abundance ra- 
tios, such as [a/Fe]. The combination of age-sensitive 
(e.g. Balmer-line indices) and metallicity-sensitive indices 
((Fe), Mgb, [MgFe]', etc.) allows in principle to derive 
accurate spectroscopic ages and metallicities. This ap- 
proach is less affected by the age-metallicity degener- 
acy than broadband photometry. In addition, combin- 
ing indices sensitive to the abundance of a-elements 
and iron ca n provide important clues on star-formation 
time scales llTinsleYlll979t lMatteuccilll993 lGreggiollT997l 
IThomas et al.lll999|) . 

The present paper makes use of the spectroscopic 
data presented in Puzia et al. (2004, hereafter Paper I) 



A fitting function gives the strength of a specific line index 
as a function of log g,T c g, and [Fe/H]. 

2 A response function gives the fractional change of an index 
as the result of an abundance change of a given element or the 
total metallicity. 
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to derive global ages, metallicities, and [a/Fe] ratios for 
globular clusters in early-type galaxies. A more detailed 
study of these global parameters as a function of internal 
galaxy properties, such as morphological type, environ- 
ment, galaxy mass, etc. will be presented in subsequent 
papers of this series. 

We describe the selection of high-quality spectroscopic 
data in The best combination of Lick indices to achieve 
most reliable age and metallicity estimates is discussed in 
f)3l In 2] we derive ages and metallicities of globular clus- 
ters in the studied early-type galaxies. Global [a/Fe] ratios 
and correlations with age and metallicity for these globu- 
lar clusters are presented in iJSJ We discuss our results in 
ffl 

2. Selection of Data 

Our entire sample contains 143 globular clusters in seven 
early- type galaxies. About half of our spectra satisfy the 
ideal S/N standards (~30 per A) to derive accurate ages, 
metallicities, and [a/Fe] ratios. It is possible to achieve an 
age resolution of ~ 1 Gyr (at ages ~ 15 Gyr) only for the 
brightest globular clusters in our sample. The typical sep- 
aration between the 15 and 14 Gyr isochrone in current 
SSP models is of the order AH/3 « 0.05 A, and - 0.1 A 
for the higher-order Balmer indices. This separation in- 
creases towards younger ages. Hence, the final sample has 
to be built from a compromise between age/metallicity 
resolution and sample size. We, therefore, set the selection 
to clusters with a statistical measurement uncertainty of 
< 0.4 A and < 0.6 A for H/3 and higher-order Balmer 
line indices, respectively. This selection corresponds to a 
minimum age resolution At/t ~ 0.3. As a metallicity in- 
dicator, we use the composite [MgFe]' index (see below). 
An error cut at 0.2 A for this index guarantees a metallic- 
ity resolution of ~ 0.25 dex at high and ~ 0.4 dex at low 
metallicities. The above selection criteria leave 71 glob- 
ular cluster spectra in our sample which correspond to 
^50% of the initial data. Henceforth, we consider only 
these selected globular clusters and combine the data of 
all sample galaxies to increase statistical significance in 
the discussed relations. Note, that this sample is biased 
towards the bright end of the globular cluster luminosity 
function, i.e. most massive clusters. Otherwise, it spans a 
wide range in globular cluster color and host galaxy prop- 
erties (see Tab.rj). 

Our original colour selection criteria prevented us from 
observing globular clusters with a combination of low 
metallicities [Z/H] < — 1.3 and relatively young ages (t SS 5 
Gyr). However, these objects are essentially non-existent 
in our sample galaxies (see Paper I). Figure ^ illustrates 
the age and metallicity ranges imposed by our colour se- 
lection. All Milky Way globular clusters would be selected 
by our colour cuts 0.8 < V - I < 1.3, 1.5 < B - I < 2.5, 
and 1.0 B — R < 1.7 (we refer to Paper I for details). 
The applied colour cuts also suggest that metal-rich old 
globular clusters might be missing in our final sample. We 
point out that, within our luminosity cut (V 23 mag), 




age [Gyr] [Fe/H] 

Fig. 1. Illustrated are colours which were used to se- 
lect globular cluster candidates as a function of age and 
metall icity and taken from model predictions of lMarastonl 
( 2005) . Upper and lower colour cuts are indicated as hor- 
izontal lines. Curves in left panels are parametrized by 
metallicity, curves in right panels are parametrized by age 
and indicated accordingly. 

the fraction of such objects amounts only a few percent. In 
order to fill the slit-masks we included cluster candidates 
redder than our colour cuts. However, only two genuine 
globular clusters were found among ~ 30 " mask- filler" ob- 
jects. Unlike the expected old ages, we find intermediate 
ages (5-10 Gyr) for these two globular clusters. The frac- 
tion of old to intermediate-age, very metal-rich globular 
clusters is likely to increase in a sample reaching fainter 
magnitudes. However, given our luminosity selection, the 
number of missed objects is negligible. 

Our dataset does not guarantee a strict consistency 
of the derived mean ages, metallicities and [a/Fe] ratios 
among the galaxies of our sample. These parameters are 
subject to change, since neither the sampled fraction of 
each globular cluster system, nor the sampling of colour 
distributions and luminosity functions is identical from 
galaxy to galaxy. Given the still relatively small num- 
bers of globular clusters per galaxy, a peculiar age and/or 
metallicity distribution in one galaxy can influence the re- 
lation between age, metallicity and [a/Fe] of the complete 
sample, although we verified that no galaxy represents a 
clear outlier (see future papers in these series). 

3. Reducing Systematic Uncertainties 

In the following we determine the best combination of in- 
dices as diagnostics for age and metallicity. Taking into 
account the uncertainty of our line index measurements, 
the mean uncertainties of the Lick system, and the limits 
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Table 1. Basic information o n ho st galaxies. A ccording to the column the references are: f 11 |de Vaucouleur s et al. 
l|l99l|, : (2) iTonrv et alJ l|200l|) : (3) iTullvl l|l988|) : (4) and (5) iKissler-Patig et al.l l|l997| ,. lAshman fc ZepJ l )l998|) . 



Galaxy 


type 


(m - M) v 


Mb 


Nac 


qb 




(1) 


(2) 


(3) 


(4) 


(5) 


NGC 1380 


-2/LA 


31.23±0.18 


-20.04 


560 ± 30 


1.5 ±0.5 


NGC 2434 


-5/E0+ 


31.67±0.29 


-19.48 






NGC 3115 


-3/L- 


29.93±0.09 


-19.19 


520 ± 120 


1.6 ±0.4 


NGC 3379 


-5/E1 


30.12±0.11 


-19.39 


300 ± 160 


1.2 ±0.6 


NGC 3585 


-5/B6 


31.51±0.18 


-20.93 






NGC 5846 


-5/E0 


31.98±0.20 


-21.16 


2200 ± 1300 


3.5 ±2.1 


NGC 7192 


-4.3/E+ 


32.89±0.32 


-20.55 







a Total number of globular clusters 

b Specific frequency, S N = N GC ■ l0 OA < Mv+1 ^ jHarris fc van den Bergblll98lH 



on the predicting power of SSP models, we construct the 
relatively best combination of diagnostic diagrams from 
Lick line indices. The found combination maximally re- 
duces the internal uncertainties of the Lick system and 
the age-metallicity degeneracy of line indices as well as 
the influence of abundance ratio variations. 

3.1. Influence of the Blue Horizontal Branch at High 
Metallicities 

As shown bv|Maraston fc Thomasl l)2000(l for H/3 and by 
iMaraston et alJ (|2003[) " for the higher-order Balmer lines, 
the morphology of the Horizontal Branch (HB), when 
extended to warm temperatures (~ 10000 K), plays a 
major role in increasing the strength of Balmer indices. 
This effect is due to metallicity and confuses the use of 
Balmer lines as pure age indicators (on t his topic see als o 
Ide Freitas Pacheco fc Barbuvlll995l and iLee et al"ll2000|) . 
Since the HB morphology cannot be predicted by first 
principles of stellar evolution, as it is determined by mass- 
loss, the line indices need to be calibrated with glob- 
ular clusters for which th e HB morphology is known 
<|Maraston fc Thomasl [20001: IMaraston et al.N2003|) . This 
exercise is clearly impossible for extragalactic globular 
clusters an d is the reason why in t he following we use both 
models bv IMaraston et afl |2003|) which include red and 
blue HB morphologies as a function of metallicity. These 
models encompass the observed range of Balmer lines 
in Milky Way globula r clust ers. As thoroughly explored 
bv iGreggio fc Renzinil ljl990|) . blue HBs are in principle 
possible also in metal-rich (Z <; Zq) stellar populations 
(see afso lRich et al.lll997t for examples of two metal- rich 
Milky Way globular clusters), that suffer from enhanced 
mass-loss and/or have a high Helium abundance. Models 
with blue HBs a t high metallicity are published elsewhere 
l|Marastonll2005|) . and here we report on their differential 
effect. 

To empirically estimate the effect of a varying hor- 
izontal branch morphology on Balmer indices at high 
metallicities, which might be responsible for the in- 



1 0.5 -0.5 -1 




1 0.5 -0.5 -1 1 0.5 -0.5 -1 

blue HBR red blue HBR red 



Fig. 2. The behavior of Lick Balmer indices as a function 
the horizont al branch morph ology (HBR) . This parameter 
is defined in lLee et all (|l994|) as HBR (B-R)/(B±V±R), 
where B and R are the number of stars blue-wards and 
red wards of the instability strip. V is the number of vari- 
able stars inside the instability strip. HBR = 1 indicates 
an entirely blue and HBR = — 1 and entirely red hori- 
zontal branch. Here we plot data f or globular cluster s in 
the Milky Way (sq uares: data from IPuzia et al.ll2002l and 
Trager ^ al.| 11998 ^ M 31 (i nverted triangles: | Puzia et alJ 
20051 iRich et alJl2003l and iTrager et al.lll998l). and th e 
Large Magellanic Cloud triangles: iBeaslev et al.ll2002b|) . 
Filled symbols show globular clusters with a metallicity 
[Z/H] -0.6 (note, all have HBR =-1). NGC 6388 and 
NGC 6441 are indicated by stars at HBR = -0.7. A small 
arrow at each panel's right ordinate indicates the most 
extreme Balmer index value for globular clusters with 
[Z/H] > —0.6. Several clusters show a tendency for higher 
Balmer indices at HBR w —0.5. Their deviations from the 
overall sequence are < 2a. 
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creased Balmer indices of metal-rich globular clusters, 
we para metrize the HB morphology with the HBR pa- 
rameter l|Lee et al.lll994|) . Figure |3| shows the strength 
of Balmer line indices as a function of HB morphol- 
ogy (HBR). In general, blue horizontal branches produce 
significantly stronger Ba lmer indices for g l obular clus- 
ters in the Milky Way llTrager et all Il998t IPuzia et all 
120021) . M3 1 3 llTrager et alJll99St IPuzia et al.ll2005|) . and 
the LMC l|Beaslev et al.ll2002bh . Each panel shows a se- 
quence of metallicity, where the scatter can be attributed 
to the "second parameter". We use the globular clus- 
ters NGC 6388 and NGC 6441 (which host the bluest 
horizontal branches among metal-rich Galactic globular 
clusters, indicated by stars in Fig. |3J) and clusters at 
[Z/H] —0.6 (filled symbols) which have entirely red HBs 
(e.g. NGC 6356 and NGC 6637) to derive a representa- 
tive "second-parameter" variation of Balmer line indices 
at high metallicities. As this approach is fully empirical 
and based on the largest HB morphology fluctuation lo- 
cally observed, we cannot rule out that even more extreme 
HB morphologies for globular clusters at a given metallic- 
ity exist outside the Local Group. We find, in the extreme 
case, offsets of 0.4 A in R(3, 3.3 A in B.'Ja, 1.4 A in H~/ F , 
2.0 A in HSa, and 1.0 A in H<5p between metal-rich glob- 
ular cluster with the weakest Balmer index (with entirely 
red HBs) and NGC 6388 and 6441. The HB morphology 
has negligible effect on the [MgFe]' index. Consequently, 
increased Balmer indices at high metallicities might well 
be, at least partly, the result of HB morphology variations. 

In the following we use the Balmer indices as age indi- 
cators with confidence at low metallicity ([Z/H] ^ —0.6), 
because HB morphology is included in our SSP models 
and under control. At high metallicity ([Z/H] —0.6), we 
have the warning in mind that ages could be degenerate 
with the presence of unresolved blue HBs. For this reason 
we will refer to such young ages as "formal" . 

3.2. The Influence of [a/Fe] Variations on Iso-Age 
Tracks 



Recently, iThomas etahl l|2003l 12004 calculated new the- 
oretical Lick index predictions which are parametrized for 
well-defined [a/Fe] ratios for a wide range of ages and 
metallicities. These models take into account the effects of 
changing element abundance ratios on Lick indices, hence 
give Lick indices not only as a function of age and metal- 
licity, but also as a function of the [a/Fe] ratio. They are 
based on the ev olutionary population synthesis code of 
iMarastonl l|l998l) . The impact fro m element-ratio changes 
is c omputed wi t h the help of the iTrinicco B"eTil l|l99fi|) 
and lKorn etaP ^05) response fun ctions, using an exten - 
sion of the method introduced by llrager et al ] fcOOOal) . 
Because of the inclusion of element-ratio effects, the influ- 
ence of [a/Fe] on Balmer indices can be studied, and is il- 



3 We derive HBR parameters for globular clusters in M31 
from colour-magnitude diagrams using HST data kindly pro- 
vided by Michael Rich (GO:6671). 



lustrated in Figure|2l Here we show the influence of [a/Fe] 
variations on frequently used age/metallicity diagnostic 
diagrams. In general, variations of Balmer-line indices for 
isochrones with [a/Fe] ratios between solar and +0.5 dex 
are of the order - 0.007-0.25 A for low, and ~ 0.18-3.63 
A for high metallicities, and introduce a relative age un- 
certainty in the range At/t ~ 0.1 — 0.2. These variations 
are due to different contaminations by metal absorption 
features inside the Balmer index passband definitions, as 
illustrated in Figure The Figure shows the higher-order 
Balmer indices generally include more metal absorption 
lines in their feature and pseudo-continuum passbands 
than the H/3 index, which is the Balmer index with nar- 
rowest passband definitions. Hence, among all Lick Balmer 
line indices, the least [q/Fel-sensitive index is H/3, fo l lowed 
byEPyp, H 7A , H(5 A , H5 F (see also IThomas et alJ l2003t 
2004). This exercise demonstrates that at high metallici- 
ties the impact of [a/Fe] variations on age/metallicity de- 
terminations (in particular, those of early-type galaxies) 
can significantly alter the results and/or introduce spuri- 
ous correlations. 



3.3. The relatively best Age Indicator 

The Balmer line series provides the best spectroscopic 
age indicator among the set of Lick line indices. The 
Lick system defines five indices ( H/3, H7A, HS\. Htf , 
and H^f) for three Balme r lines jWorthev et alJ [l994; 
IWorthev k. Ottavianilll997|) . Figure 0] shows the passband 
definitions for all Balmer indices. In combination with a 
metallicity diagnostic, these higher-order Balmer line in- 
dices are widely used to determine (luminos ity- weighted) 
ages and me t allicities of galax i es (e.g. | Tra ger et alJ 
I l998t l2000albl iKuntschnerl l200(i IPoggianti et all l200li 
Kuntschner et alJl2002al IThomas et alJl2003|) . ~ 

However, different types of diagnostic plots based on 
different Balmer-line indices are used throughout the liter- 
ature. Although the age predicting power of an arbitrarily 
chosen diagnostic plot (most common versions include the 
H/3 and Mg2 or (Fe) indices might yield accurate-enough 
results for a specific scientific goal (e.g. the mean age dif- 
ference between two different galaxy samples), the choice 
of a specific diagnostic plot is still subject to observational 
constraints and individual assessment, and makes com- 
parisons between studies difficult. As a consequence, most 
authors use several diagnostic plots with different Balmer 
line indices and assign equal importance to the results de- 
rived from each of those. 

In the following we provide a recipe to define a quantity 
from which the relatively best Balmer-line age indicator 
can be determined. This quantity takes into account the 
quality of a given dataset and the diagnostic power of the- 
oretical predictions from which one intends to derive the 
age and metallicity. It does not take into account system- 
atic uncertainties in the fitting functions for a particular 
index. 
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Fig. 3. Infl uence of [a/Fe] variations on age-metallicity diagnostic grids. Panel (a): H/3 vs. (Fe). Lines show model 
predictions l|Thomas et alJEoO.ll 12004^ for stellar populations with metallicities [Z/H]=— 2.25 to 0.67 dex and two ages 
13 (solid lines) and 3 Gyr (dot-dashed lines) parametrized for three different [a/Fe] ratios 0.0, 0.3, and 0.5 dex. Panel 
(b) shows the same model predictions for H/3 vs. [MgFe]', while panel (c) illustrates the changes in the H/3 vs. Mgb 
grid. The other panels show the predicted variations for higher-order Balmcr line agc/mctallicity diagnostic grids. 



In particular, the age sensitivity of an index is a func- 
tion of the following parameters: 

• 7/: mean error of the data in A 

• transformation accuracy to the Lick system in A 

• 7: mean error of the original Lick spectra in A 

• 5: accuracy of the Lick fitting functions in 
AllWorthev et alJlT^ IWorthev fc Ottavianilll 997h 

• T>z'- index range in A covering all ages at a given metal- 
licity, hereafter termed the dynamic range 

• Sa,z,t'- degeneracy parameter, which quantifies the sen- 
sitivity of an index / to age, metallicity, and a/Fe ratio 
at a given [a/Fe], metallicity, and age (i.e. the impact 
of the age-metallicity and age-a/Fe degeneracy), given 
in dex/ Gyr 



The numerical values for each parameter are given in 
Table [21 for each Balmer index. It is worth noting that 
some of these values are only valid for our data qual- 
ity in combi nation with the SSP models of lThomas et alJ 
1I2OO3U2OO4I) . For different data and SSP models, 77, V z , 
and <So,,z,t are subject to change. To quantify the most 



age-sensitive and least metallicity-sensitive Balmer index, 
we define the quantity 



n 



• <?a!,Z,i 



vV + C 2 + 7 2 + $ 2 
where the degeneracy parameter, <S Q ,z,t, is defined as 



(1) 



dl 




dl 


2 

+ 


dl 


') 


~dt 




<9[Z/H] 


a,t 


d [a/Fe] 


z,t / 



(2) 



1Z, which is given in dex- Gyr -1 , is essentially the dy- 
namic range of an index / at a given age and metallicity 
expressed in units of the total uncertainty. We use the 
mean dynamic age range 2?z at two different metallicities 
[Z/H]= —1.35 and 0.0 between the 1 and 15 Gyr isochrone. 
Each SSP model provides a well-defined relative age scale 
which we use here to parametrize iS Qi z,t for two different 
metallicities [Z/H]= —1.35 and 0.0, two different ages 3 
and 13 Gyr, and [a/Fe] = 0.3. 5 a ,z,t is the ratio of age, 
metallicity, and [a/Fe] partial derivatives (see Eqn.[5J. In 
other words, 5 Qi z,t is a measure of the age-metallicity and 
age-a/Fe degeneracy, and is maximal for indices which are 
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Fig. 4. Passband definitions for Balmer-line, Mgb, Fe5270, and Fe5335 Lick indices with their feature and adjacent 
continuum passba nds. The over-plotted spectrum is a high-S/N spectrum of the Galactic globular cluster NGC 6284 
llPuzia et alJl2002l). The resolution is ~ 7 A and was left untouched to keep satellite lines visible. Line data were taken 
from iReader fc Corliss! (ll98]J) . data for molecular bands are from iRosenl l|l952t) . Note the large number of satellite 
lines which are included in the passband definitions. 



most sensitive to age and least sensitive to metallicity and 
[a/Fe] variations, at the same time. 

The highest 1Z indicates the best age indicator with 
least age-metallicity and age-a/Fe degeneracy. In Table 
we provide values for T>z at two different metallicities and 
for <5>Q,,z,t at [a/Fe] = 0.3 and four age-metallicity combi- 
nations for each Balmer line index. Since SSP models do 
not provide continuous but discrete predictions the par- 
tial derivatives are substituted by difference ratios, e.g. 
dl/dt — > AI / At. The quotients are determined by linear 
interpolation of SSP models. 

We determine the relatively best age indicator from 
the set of five Lick Balmer indices by combining the mean 
dynamic range (2?z), the mean age-metallicity and age- 
a/Fe degeneracy parameter (<S Qi z,t), and the total index 
uncertainty which is the denominator in equation ^ The 
final mean 1Z is documented in the last column of Table [5] 
We find that the relatively best age diagnostic for our data 
is the H7A index followed by the indices H/3 and H6a- Hyp 



and HSp have the smallest 1Z values and are not considered 
to be reliable age indicators. 

It is instructive to see that despite the relatively large 
age-metallicity degeneracy of the H7A index, it still pro- 
vides the most accurate age predictions. This fact is pri- 
marily due to the large dynamic range of H7A compared 
to its mean measurement uncertainty. H/3, on the other 
hand, has a relatively large total uncertainty and the mea- 
surements will therefore be more scattered over the di- 
agnostic plot's parameter range. However, H/3 shows by 
far the highest {S a> z,t) value (column 13 in Tab. [2 see 
also Fig. [3J) and is therefore least sensitive to the cumu- 
lative effects of age-metallicity and age-a/Fe degeneracy. 
In general, the higher-order Balmer lines require less S /N 
to guarantee a similar total index accuracy as H/3, which 
is primarily due to the narrow passband definition of H/3 
(see Fig. QJ. If our dataset would be infinitely accurate 
(i.e. 1) = in Eqn. QJ, the order of 1Z from the best to 
worst Balmer index would remain unchanged. The value 
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of 1Z is predominantly governed by uncertainties in the 
fitting functions of the respective index. To vary this or- 
der, the mean measurement uncertainties have to be very 
discrepant and the SSP-modcl predictions have to devi- 
ate significantly from the model used here. It is expected 
that the relative accuracy of Balmer index measurements 
is comparable between different datasets as they are usu- 
ally derived from one optical spectrum. The relative age 
scale of SSP models appears to be quite stable against the 
choic e of different stellar e volutionary tracks for ages > 1 
Gvr JCharlot et al.lll99ri iTrager et all l2000al: iMarastonl 
120031) . This scale is used in our above prescription. It can 
therefore be expected that no large fluctuation in 1Z will 
arise from the use of different SSP model predictions 4 . 

3.4. The relatively best Metallicity Indicator 

The index with the highest metallicity sensitivity and min- 
imal age-sensitivity could in principle be found in a com- 
parable way as it was done for the relatively best age di- 
agnostic. The major impact of typical metallicity tracers, 
such as (Fe), Mg2, and Mg6, on the absolute metallicity 
scale is expected to arise from changing abundance ratios. 
To reduce the influence of Ja/Fe] variations on age and 
metallicity determinations, iThomas et alJ l)2003|) modify 
the [MgFe] index 5 to obtain an [a/Fej-insensitive metal- 
licity index, 



[MgFe]' = x/Mgb ■ (0.72 x Fe5270 + 0.28 x Fe5335). 

Figure 01 shows the behavior of isochrones in three dif- 
ferent frequently used diagnostic plots. The middle column 
impressively illustrates that [MgFe]', indeed, is essentially 
independent of [a/Fe]. Henceforth, we adopt the [MgFe]' 
index as the best metallicity indicator and use it in combi- 
nation with selected age indicators (see Sect. 14. l|) to derive 
ages, metallicities, and [a/Fe] ratios for our sample glob- 
ular clusters. 



3.5. Excluding contamination by Ionized Gas 

It is known that ~ 40 — 60% of early- type galaxie s show in- 
dications of emission in their absorption spectra f Caldwel 



dications ot emission m their absorption spectra (L/aldwel 
1984 [IPhillips et al.ll98fitfGonzalezll99ilrxOudfrooii et al 



1994). In a narrow-band imaging survey. iMacchetto et al 
(1996J) find ionized gas in ~ 80% of early-type galax- 
ies including flocculent Ha+[Nn] emission in NGC 3379, 
NGC 5846, and NGC 7192, well within - li? eff , but 
no significant emission in NGC 3115. This gas is located 
in the central parts and distributed in a rather regular 
way, suggestive of a disk. If dominant, all Balmer in- 
dices would be affected, along with potential contami- 
nation of Fe5015 by [Qui] (A 5007 A] and o f Mg6 by 
[Ni (A 5199 A) llGoudfrooii fc Emsellemlll 996ft . In that 

4 We note that systematic uncertainties in the fitting func- 
tions and effects of emission-line filling are not considered by 
this exercise. 

5 [MgFe] = y/Mgb- (Fe), see lOonzalej l)l99.^ . 



case, our measurements would indicate too old ages. This 
effect rapidly dec reases from H/3 towards H7 and US 
l|Osterbrookl lT989). i.e. higher-order Balmer indices are 
less affected and should be preferentially used for age de- 
terminations in the presence of ionized gas. 

In order to exclude a major effect of ionized gas on 
line-strength measurements in our globular cluster data, 
we performed several tests/estimates. 

Since lin e emission is concentra ted in the central parts 
of galaxies l|Macchetto et al.lll996h . we expect a correla- 
tion of Balmer indices with galacto-centric radius if line- 
emission contamination is significant. We find no evidence 
that Balmer indices are correlated with galacto-centric 
distance. A more detailed analysis of background spectra 
shows that most clusters are located within ^2 — 3 R c g 
and that the flux level of the diffuse galaxy light is well 
below the object flux. In particular, we find no correlation 
inside one effective radius, where line emission is expected 
to be strongest. Furthermore, we find no correlations of 
Balmer indices for globular clusters with Balmer indices 
measured on corresponding background spectra. 

Visual re-inspection of the background subtraction 
process for some low-H/3 outliers (see Fig. EJ underline 
the good quality of background modeling and subtraction. 
However, problems with accurate background subtraction 
might occur in very few cases when line emission has a 
very filamentary structure. For instance, the worst case 
scenario would be when a globular cluster overlaps with a 
filament of ionized gas while the slit is aligned perpendic- 
ular to such a filament. H owever, a filamentary em ission 
pattern is not found in the IMacchetto et alJ l)l996|) study 
for our host galaxies. We conclude that line emission has 
no measurable effect on the scatter in the age/metallicity 
and [a/Fe] diagnostic plots. 



4. Ages and Metallicities 

4.1. Iterative Method 

Having determined the relatively best metallicity diagnos- 
tic, we determine the best combination of Balmer-line in- 
dices as our prime age indicator. In deciding whether a 
Balmer index will be chosen as part of this most reliable 
age proxy, we carefully analyze the following points for 
each Balmer index and assign priorities in descending or- 
der: 

a) ranking of the reliability parameter 1Z, as described in 
Section 13. 31 

b) including the maximum number of the Balmer lines 
without duplicate measurements by two or more line 
indices 

We choose the combination of H7A, H/3, and H(5a as our 
most reliable age indicator. Consequently, ages and metal- 
licities for individual globular clusters are computed as the 
weighted mean of the parameters derived from diagnos- 
tic plots constructed from these three Balmer indices vs. 
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Table 2. Summary of the coefficients relevant to equation ^ The coefficients in columns 2 — 8 are given in units of 
A. Columns 9 — 14 are given in dex/Gyr. Note that for the higher-order Balmer indices column 5 giv es the arithmetic 
mean of the uncertainties for cool and warm stars (see Table 3 and 4 in lWorthev fc Ottavianiiri997|) . 
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[MgFe]'. As weights we choose the reliability parameters 
Tl from Table 

Since the SSP model predictions allow full control over 
[a/Fe] variations within the diagnostic grids, we use an it- 
erative approach in combination with a second diagnostic 
grid from which we derive the [a/Fe] ratios (Mg2 vs. (Fe), 
see Sect. |3J. As a first step the a-enhancement for each 
individual globular cluster is derived. This value is used 
to interpolate the age/metallicity diagnostic grid for the 
correct a-element enhancement. From so adjusted grids 
ages and metallicities are computed using linear 6 interpo- 
lation of the model predictions employing a least-square 
technique. These ages and metallicities are then used to 
adjust the [a/Fe] diagnostic grid, which itself is slightly de- 
pendent on age. This procedure is iterated until the age, 
metallicity, and [a/Fe] values converge. Extensive tests 
show that this goal is reached within a few iteration steps. 
Along with these values we also determine the statistical 
errors due to index measurement uncertainties. Typically 
the errors are asymmetric which is a result of the skewness 
of diagnostic grids (see Fig. EJ . 

For objects which fall outside the diagnostic grid, we 
assign the most likely extreme grid value and do not ex- 
trapolate the model predictions, i.e. we project the data 
onto the closest grid point along its error vector 7 . Note 
that due to the increasing influence of hot blue horizontal- 
branch stars, the iso-age tracks for old (t > 8 Gyr) stel- 
lar populations with metallicities below —1.0 dex tend to 
overlap. This introduces an ambiguity in assigning ages 
and metallicities to individual globular clusters and artifi- 
cially broadens the age distributions at old ages. However, 
only six objects are affected by this ambiguity which has 
no effect on the following results. For these six objects the 
routine randomly assigns either young or old ages and the 
corresponding metallicities. 

4.2. Consistency Checks 

The comparison of derived formal ages and metallicities of 
individual globular clusters shows that different diagnos- 
tic diagrams make age predictions which are somewhat 

6 Note that the SSP model grids are interpolated in linearly. 

7 To illustrate the effect of this procedure we mark objects, 
which fall off the grid in all diagnostic plots, by double-hatched 
histogram in the distributions plots in Figure |7] 



inconsistent with one another. The average relative age 
uncertainty varies in the range At/t w 0.2 — 0.5 which we 
attribute to systematic errors in the age scale of model 
predictions and/or systematics in the calibration of the 
data. By combining age predictions from different diag- 
nostic plots we smooth out these systematic effects which 
influence age determinations derived from single Balmer 
indices. 

This age inconsistency is partly surprising as the 
Maraston et al. (2003) models are well calibrated for Milky 
Way globular clusters and their ages derived from different 
Balmer indices are consistent with each other. 

On the other hand, we find very good agreement of 
metallicity predictions from all diagnostic plots within an 
average uncertainty of ~ 0.2 dex. 

Note also, that a striking feature of all age/metallicity 
diagnostic plots in Figure [S] is that a significant fraction 
of globular clusters lie below the oldest iso-age track. This 
is observed in other samples of e xtraqalactic globular clus- 
ters in the literature as well (e.g.lKissler-Patig et al.lll998t 
ICohen et al]ll99St IPuzia et al.ll2000t lLarsen et al.ll2003|) 
Systematics errors in the data reduction have been deemed 
extremely unlikely given the good Lick standard calibra- 
tion (see Paper I) and are also unlikely given the general 
problem in the literature. We could also exclude emission 
filling as source for the Balmer index inconsistencies which 
is presented in Section l3~5l 

To test whether the spectroscopic age and metallicity 
predictions for these outliers are consistent with their ob- 
served photom etric colors (see P aper I for details) , we use 
SSP models of iMarastonl l)2005|) to compute photometric 
colors from the derived spectroscopic ages and metallici- 
ties. Figure shows the comparison of computed and ob- 
served data. We find good agreement for the B— V, V—I, 
and V — K colors. We find no systematics in color resid- 
uals as a function of Balmer index with respect to the 
one-to-one relation. However, we find significantly smaller 
residuals towards more metal-rich objects. 

We also investigate if higher-quality data alleviates 
these inconsistencies and select a sub-sample with sta- 
tistical errors AH/3 < 0.11 A and AH 7i<5 < 0.165 A for 
higher-order Balmer lines, which is shown as solid dots 
in Figure |5] For this sub-sample, we indeed find a better 
agreement between computed and observed photometric 
colors. We take this as an indication that the age and 
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Fig. 5. Computed vs. measured photometric colors for our sample globular cl u sters. Computed colors were derived 
from spectroscopic ages and metallicities using the SSP models of iMarastonl l|2005(l . The circles show all selected 
globular clusters for which optical/near-infrared colors were obtained (see Paper I for details). Solid circles show the 
high-quality sample, which was selected by a more constrained error cut (see text for details). Errorbars indicate 1-a 
uncertainties. 



metallicity predictions for the higher-quality dataset are 
more robust. Based on our photometric cross-check, we 
only consider the ages and metallicities derived for our 
high-quality sub-sample as trustworthy. However, to illus- 
trate the difference between the high-quality and the rest 
of the data we plot the remaining dataset in the following 
diagnostic plots as open symbols in Figure and show 
their distributions as hatched histograms in Figured 

4.3. Results 

In the following, we discuss ages and metallicities for 
the high-quality sub-sample derived from diagnostic plots 
using the Balmer-line indices H/3, H7A, and US a (see 
Sect 14. 211 . The corresponding diagnostic diagrams are 
shown in Figure along with the other two diagrams con- 
structed from the Balmer indices EPyp, and HSp. Since for 
seven globular clusters one or more of the three Balmer 
indices H/3, H7A, and H<5a are not available, we exclude 
these objects from the subsequent analysis to avoid sys- 
tematics, which leaves 17 globular clusters in our high- 
quality sub-sample. Their age and metallicity distribu- 
tions are shown in Figure as solid histograms. 

Our method reveals that a significant fraction of clus- 
ters in the high-quality sample (6/17 or 35%, see solid his- 
togram in Fig. 0) has formal ages between 5 and 10 Gyr. 
Only one object has a formal age below 5 Gyr (~ 3.8 Gyr) 
and there are no globular clusters with extremely young 
ages (^ 1 Gyr). We split the sample at [Z/H] = —0.6, 
corresponding to the dip in the Milky W ay globular clus- 
ter metallicity distribution l)Harrislll996|) . into metal-poor 
and metal-rich globular clusters. With this distinction we 
find an increase in age spread from metal-poor to metal- 
rich clusters. We derive a dispersion of 0.5 Gyr for metal- 
poor and a dispersion of 3.5 Gyr for metal-rich globular 
clusters. A weighted linear least-square fit to our high- 
quality dataset reveals a weak age-metallicity relation, in 
the sense that more metal-rich globular clusters appear on 



average younger. The significance of this relation needs to 
be tested with larger samples. 

Our high-quality globular cluster sample covers metal- 
licities in the range —1.3 [Z/H] ^ +0.5, with a mean 
—0.36 ± 0.13 and dispersion a = 0.52 dcx. Four out of 17 
globular clusters have formally super-solar metallicities. 
All clusters that are younger than 10 Gyr have metallic- 
ities [Z/H] <; —0.4, which excludes their formation from 
primordial gas clouds. 

We note that no super metal-rich, intermediate-age 
counterpart glob ular cluster p opulation is found in the 
Milky Way (e.g. Harris 2001). However, globular clus- 
ters with intermed i ate ages are found in M31 (e.g. 
Barmbv et all l200(j iBurstein et all 120041: iBeaslev et all 
2004U P uzia et al 2005), whose bulge-to-disk ratio is larger 
than that of the Milky Way. An interesting idea is that 
this mode of globular cluster formation may not have man- 
ifested iteelf_2n_oui_G^laxvdue to Milky Way's smaller 
bulge l|Goudfrooii et alJl2003j) . 

In summary, the conclusions that can be drawn from 
our previous analysis of globular clusters in early-type 
galaxies are the following: 

• metal-poor globular clusters are on average older than 
metal-rich globular clusters 

• there is a tendency for metal-rich globular clusters to 
have a younger mean age relative to their metal-poor 
counterparts 

• metal-rich globular clusters exhibit on average a larger 
age spread than metal-poor globular clusters 

5. [a/Fe] Ratios 

In this Section we derive [a/Fe] ratios for our high-quality 
sample globular clusters using a diagnostic diagram which 
is least sensitive to age/metallicity variations. Such a dia- 
gram can be constructed from the indices (Fe) and Mg2, 
which primarily trac e the abundances of i ron and the a- 
element magnesium llTrinicco Belll Ffl95) . We note that 
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Fig. 6. Age-metallicity diagnostic plots (panels a~e) constructed from Balmer indices H7A, H/3, HS\, fPyp, and 
H5f vs. [MgFe]' for high-quality glob ular cluster spectra (solid dots). Open circles show our remaining selected 
globular cluster data. SSP models of iThomas et alJ (|2004f) have been over-plotted for [a/Fe] = 0.3, metallicities 
[Z/H]=-2.25, -1.35, -0.80, -0.33, 0.00, 0.35, and 0.67 dex (dotted lines), and ages 15, 10, 5, 3, 2, 1 (solid lines). 
The red dotted line is an interpolated iso-metallicity track for [Z/H] = —0.8 and is used to split the sample between 
metal-poor and metal-rich globular clusters. The blue solid line is the 5 Gyr isochrone, and is used to split between 
old and young globular clusters. A dashed isochrone illustrates the Balmer index strengths for a 15 Gyr old stellar 
population with an entirely red horizontal-branch morphology. The [a/Fe] d iagnostic plot (panel f) shows the mean 
iron index (Fe) as a function of Mg2. SSP models of IThomas et al.l (2003, 200J) with constant [a/Fe] ratios have been 
over-plotted for [Z/H] between —2.25 and +0.67 dex and two ages 13 (solid lines) and 3 Gyr (dotted lines) with various 
[a/Fe] ratios 0.0, +0.3, and +0.5 dex. 
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Fig. 7. Histograms for ages, metallicities, and [a/Fe] ratios for globular clusters in early- type galaxies. The panels in 
the left column show distribution of each parameter for all selected (hatched histograms), outliers in diagnostic plots 
(double- hatched histograms), and high-quality data (solid histograms). Each panel in the right column illustrates a 
correlation of two of these parameters for the high-quality sub-sample. Dashed lines in the correlation plots show 
weighted least-square fits. 



among the three Mg-sensitive indices, Mgi, Mg 2 , and Mg6, 
theoretical index predictions show a relatively large spread 
in (Fe) and Mg2 for [a/Fe] ratios between solar and ~ 0.5 
dex, at high mean metallicities. Given the quality of our 
data we can expect a good discrimination between en- 
hanced and solar- type [a/Fe] ratio for individual globular 
clusters at metallicities [Z/H] <; —0.8 and a good estimate 
of the mean [a/Fe] ratio at lower metallicities. 

Figure HJ (panel /) shows that the Mg2 vs. (Fe) diag- 
nostic diagram is not entirely free from the age/metallicity 
degeneracy. Iso-[a/Fe] tracks for three different ratios (0.0, 
0.3, and 0.5 dex) are plotted for two ages (3 and 13 Gyr, 
indicated by dotted and solid lines, respectively). It is ob- 
vious that age information is needed to choose the correct 
set of tracks for a reliable [a/Fe] determination. The diag- 
nostic grid is interpolated in our iterative fitting routine 
using information derived from age/metallicity diagnostic 



grids (see Sect. 14. If) . In other words, [a/Fe] ratios are de- 
termined simultaneously with ages and metallicities. The 
averaging is identical to the one for ages and metallicities 
m Section Ol 

A histogram of [a/Fe] ratios for our high-quality sam- 
ple is shown in Figure All globular clusters are consis- 
tent with super-solar [a/Fe] ratios. The mean [a/Fe] of 
the sample is 0.47 ± 0.06 dex, with a dispersion of 0.26 
dex. 

Using a weighted linear least-square fit, we find ev- 
idence for a [a/Fe]-metallicity relation in the sense that 
more metal-rich globular clusters have lower [a/Fe] ratios. 
However, due to the reduced [a/Fe] resolution of Lick in- 
dices at low metallicities and our modest sample size (es- 
pecially at the metal-poor end), this trend needs to be 
confirmed with more data. Consistent with the previously 
found age-metallicity and [a/Fe]-metallicity relation, we 
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find no evidence for a [a/Fe]-age correlation. Globular 
clusters at all ages appear to have on average super-solar 
[a/Fe] ratios. 

Recent measurements of [a/Fe] ratios in globular clus- 
ter systems in other early-type galaxies reveal very similar 
results. For instance, the data of lKuntschner et al.l |2002bl 
for globular clusters in NGC 3115 show that most clusters 
are consistent with [q /Fe] ~ 0.3 oyer the e ntire range of 
sampled metallicities. lLarsen et, all l|2002bh find a super- 
solar mean [a/Fe] ratio of +0.4 dex for globular clusters of 
all metallicities in NGC 4594 (Sombre ro). Using SSP mod- 
els with non-constant [a/Fe] ratios. iForbes et"afl l)200l[) 
argue that at least some globular clusters in NGC 1399 ex- 
hibit super-solar [a/Fe] ratios. However, the found [a/Fe]- 
metallicity correlation is seen in this sample for the first 
time. This is mainly due to the higher quality of our data 
and a more reliable age/metallicity determination com- 
pared to previous studies. 

Note that super-solar [a/Fe] ratios have also been re- 



ported for the diffuse liqht of early-type galaxies (e.g. 


Davies et al.ll993UTrager et al.l2000bll 


Davies et al.ll200ll 


Thomas et al.l 120051 Kuntschner et alJ 


20023). In subse- 



qucnt papers of this series we will compare the index mea- 
surements of globular clusters with indices derived from 
our long-slit spectra of the hosts' diffuse stellar light. 



In the following we summarize the major findings of 
this section: 

• globular clusters in early-type galaxies have on average 
super-solar [a/Fe] ratios with a mean ~ 0.45 dex 

• metal-rich globular clusters have on average lower 
[a/Fe] ratios than their metal-poor counterparts 

• young and old globular clusters have on average similar 
[a/Fe] ratios 

6. Discussion 

6.1. Assembly History of Early-Type Galaxies 

If the dispersion in Balmcr-[MgFe]' diagnostic plots is en- 
tirely driven by age 8 , we find indications for a significant 
fraction of relatively young globular clusters in early-type 
galaxies. About ~ 1/3 of our sample globular clusters have 
formal ages younger than 10 Gyr, implying formation red- 
shifts Zf ^ 1.7 (in a ACDM concordance universe). 

In Figure|H|we compare the derived ages and metallici- 
ties of our high-quality sample with the corresponding dis- 
tributions of globular cluster syste ms simulated accordin g 
to a hierarchical merging scenario (Be aslev et al.ll2002a|) . 
As the average luminosity of our sample host galaxies is 
(-MB)sampic = -20.4±0.6 (see Tab.HJ), we use the globular 
cluster system of a typical simulated galaxy with a similar 
lumino sity Mb = —20.42 (see also Fig. 6 in lBeaslev et alJ 
2002a) for comparison. We split both cluster samples (sim- 
ulated and observed) in metal-poor and metal-rich clus- 
ters at [Z/H] = -0.9. The top right panel of Figure El 

8 A population of warm horizontal branch stars remains a 
potential cause for increased Balmer line strengths. 



shows that our sample is biased towards metal-rich glob- 
ular clusters, and the comparison of metal-poor globular 
clusters is limited to the statement that both observed and 
simulated clusters are all consistently old. For metal-rich 
globular clusters there is good agreement in the age range 
covered by data and simulations. 

The somewhat naive prediction of more extended for- 
mation timescales of more massive structures in the hier- 
archical picture holds only if the fraction of gas-poor to 
gas-rich mergers, the so- called dry to mixed merger r atio, 
is constant with redshift. lKhochfar fc Burkertl l|2003lh pre- 
dict that this ratio depends on galaxy mass in the sense 
that most massive ellipticals formed early in rather dis- 
sipationl ess (dry) mergers of bulge-dom inated precursors 
(see also iKauffmann fc Haehnehi l2000|) . This would im- 
ply that on average low-luminosity ellipticals would har- 
bor a higher fraction of young globular clusters, while 
most massive galaxies would preferentially host old glob- 
ul ar cluster systems. Th is is reflected in the simulations 
of iBeaslev et all (|2002aj) . The lower panels of Figure 00 
illustrate the mean ages and metallicities of globular clus- 
ter systems in faint and bright galaxies with a cut at 
Mb = —21.0. The plots show that while host galaxy lu- 
minosity affects the mean metallicity of globular cluster 
systems only marginally, a significant fraction (~ 25%) 
of low-luminosity galaxies is predicted to host globu- 
lar cluster systems with an average age younger than 
~ 10 Gyr, compared to a tiny fraction (~ 2%) in bright 
galaxies. Unfortunately, the transition between dry and 
mixed merger-dominated evolution is predicted to occur 
for early-type galaxies in the luminosity range between 
Mb ~ —20 and —21 mag, where all of our sample galax- 
ies reside. 

To first order, previous spectroscopic studies find evi- 
dence for mostly old globular c luster systems in the mas - 
sive Fornax galaxy NGC 1399 faissler-Patig et alJll998h . 
and t he two Virgo galax ies M87 l|Cohen et al.lll998i) and 
M49 llCohen et al.ll2003h . which appears to fit into the 
framework of dissipationless merging with an early assem- 
bly. We will get back to this point in a subsequent paper 
in which globular cluster systems will be investigated as a 
function of the host galaxy property. 

6.2. Formation Timescales 

Stellar populations with super-solar [a/Fe] ratios, as they 
are observed in massive elliptical galaxies, are generally in- 
terpreted as the result of very short formation time scales. 
However, in the case of dissipative merging, hierarchical 
merging models predict frequent merging events, which 
is expected to resu lt in relatively low [a/Fe] ratios (e.g. 
iThomas et al . 1999). Although somewhat on the high side, 
our formal mean [a/Fe] ratio of ~ 0.47 dex for globular 
clusters in early-type galaxies is in line with values mea- 
sured for the diffuse light. Hence, the formation time scales 
of field stars and globular clusters appear to be similar in 
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Fig. 8. Top row: Comparison of globular cluster ages and metallicities in our sample galaxies with those in a simulated 
host galaxy with a similar luminosity (mass). All globular cluster in our sample are shown as hatched histograms, 
outliers in diagnostic plots are indicated by double-hatched histograms (see Sec. l4.2|> . and high-quality data is shown as 
solid histograms. Observed and simulated globular cluster samples were split into metal-poor and metal-rich globular 
clusters at [Z/H] = —0.9 and the shading is accordingly transfered to the age distribution plot. Bottom row: Average 
ages and metallicities of sim ulated globular cluster systems in faint and bright host galaxies. All predictions were taken 
from lBeaslev et alJ (|2002al) . 



elliptical galaxies and not immediately compatible with 
hierarchical merging models. 

The found evidence for a [a/Fe]-metallicity relation 
in early-type globular cluster systems is consistent with 
a generic chemical enrichment scenario in which more 
metal-rich stellar populations have lower a-element en- 
hancements. Two types of supernovae contribute to the 
build-up of a-peak (type II) and iron-peak elements (type 
la) in the interstellar medium on diffe rent time scales (e.g. 
iMatt eucci ll994llThomas et alJll99fl|) . because of the evo- 
lutio nary delay of ~ 1 Gyr of type-la supernova progeni- 
tors l)Greggiolll997f) . Since globular clusters are not mas- 
sive enough to support significant self-enrichment, their 
[a/Fe] ratios reflect the large-scale chemical conditions 
during their formation. However, to allow a detectable in- 



crease in [a/Fe], for high-metallicity stellar populations 
it requires a larger number of type II supernovae to out- 
weigh any previous metal-enrichment by type la super- 
novae compared to lower metallicities. This implies that 
high [a/Fe] ratios at high metallicities are likely the re- 
sult of shorter formation timescales and/or higher star- 
formation rates than similar [a/Fe] ratios at lower metal- 
licities. Careful modeling of the chemical evolution of a 
globular cluster system is needed to quantify these forma- 
tion timescales. The relation between [a/Fe] and metallic- 
ity found here provides the first important constraint for 
this exercise. 

It is also important to test whether a steeper [a/Fe]- 
age relation can be found for globular cluster systems in 
more massive galaxies, which would be expected if these 
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systems formed on more extended timescales, as predicted 
by hierarchical merging models. 

7. Summary and Conclusions 

We have conducted a study of ages, metallicities, and 
[a/Fe] ratios of extragalactic globular clusters in early- 
type galaxies, based on the Lick index system. We find 
that up to ~ 1/3 of our sample of globular clusters have 
ages formally younger than 10 Gyr. This result is not bi- 
ased by one single globular cluster system in our galaxy 
sample and appears representative for early-type galax- 
ies in general. We cannot state with confidence whether 
the younger ages (i.e. t < 10 Gyr) are real or due to an 
unexpected blue horizontal branch morphology at high 
metallicities. If the younger ages are real, the found frac- 
tion should be taken as an upper limit, since our data 
only probe the bright end of the globular cluster luminos- 
ity function where relatively young globular clusters, if 
present, are expected to reside. In general, less than ~ 1/3 
of the brightest ~ 10% of globular cluster systems in early- 
type galaxies which we sample with our study could have 
formed at redshifts zj 1.7. 

We find that the formal age scatter increases and the 
mean age decreases from metal-poor to metal-rich glob- 
ular clusters, resulting in an age-metallicity relation. For 
our high-quality sample globular clusters with metallici- 
ties [Z/H] < —0.6 we find an age dispersion of 0.5 Gyr. 
The metal-rich sub-sample has a dispersion of 3.5 Gyr. 

Our high-quality sample spans a wide range in metal- 
licity between [Z/H] « —1.3 and ~ +0.5 dex. We find ev- 
idence for an [a/Fe]-metallicity relation in the sense that 
more metal-rich globular clusters have lower a-element en- 
hancements, which needs to be confirmed with more data 
that better sample the metal-poor regime. However, there 
is no indication for an [a/Fe]-age correlation. 

[a/Fe] ratios are found to be on average super-solar 
with a mean 0.47 ±0.06 dex and a dispersion of ~0.3 dex, 
which indicates formation timescales shorter than ~ 1 Gyr. 
In other words, the progenitor clouds of globular clusters 
were predominantly enriched by type-II SNe. 

Acknowledgements. We are grateful to Michael Rich for pro- 
viding colour magnitude diagrams of M31 globular clusters 
prior to publication. Many thanks go to Michael Beasley for 
sending electronic tables of his hierarchical-clustering simula- 
tions. We also thank Scott Trager for a very constructive ref- 
eree report. THP gratefully acknowledges the support by the 
German Deutsche Forschungsgemeinschaft, DFG project num- 
ber Be 1091/10-2, and the support in form of an ESA Research 
Fellowship. 

References 

Arimoto, N. & Yoshii, Y. 1987, A&A, 173, 23 
Ashman, K. A., Bird, C. M., & Zcpf, S. E. 1994, AJ, 108, 
2348 

Ashman, K. M. & Zcpf, S. E. 1998, Globular Cluster 
Systems, Cambridge University Press 



Barmby, P., Huchra, J. P., Brodie, J. P., Forbes, D. A., 
Schroder, L. L., & Grillmair, C. J. 2000, AJ, 119, 727 

Baugh, C. M., Cole, S., Frenk, C. S., & Lacey, C. G. 1998, 
ApJ, 498, 504 

Beasley, M. A., Baugh, C. M., Forbes, D. A., Sharpies, 

R. M., & Frenk, C. S. 2002a, MNRAS, 333, 383 
Beasley, M. A., Hoyle, F., & Sharpies, R. M. 2002b, 

MNRAS, 336, 168 
Beasley, M. A., Brodie, J. P., Strader, J., Forbes, D. A., 

Proctor, R. N., Barmby, P., & Huchra, J. P. 2004, AJ, 

128, 1623 
Bender, R. 1988, A&A, 202, L5 

Bender, R. 1996, IAU Symp. 171: New Light on Galaxy 

Evolution, 171, 181 
Bender, R., Ziegler, B., & Bruzual, G. 1996, ApJL, 463, 

L51 

Bernstein, R. A. & McWilliam, A. 2002, IAU Symposium, 
207, 739 

Bower, R. G., Lucey, J. R., & Ellis, R. S. 1992, MNRAS, 
254, 601 

Brodie, J. P. & Huchra, J. P. 1991, ApJ, 379, 157 
Burstcin, D., Faber, S. M., Gaskell, C. M., & Krumm, N. 

1984, ApJ, 287, 586 
Burstcin, D., ct al. 2004, ApJ, 614, 158 
Caldwell, N. 1984, PASP, 96, 287 

Carretta, E., Gratton, R. G., & Sncdcn, C. 2000, A&A, 
356, 238 

Chariot, S., Worthcy, G., & Bressan, A. 1996, ApJ, 457, 
625 

Cohen, J. G., Blakeslee, J. P., & Ryzhov, A. 1998, ApJ, 
496, 808 

Cohen, J. G., Blakeslee, J. P., Cote, P. 2003, ApJ, 592, 
866 

Cole, S., Lacey, C. G., Baugh, C. M., & Frenk, C. S. 2000, 

MNRAS, 319, 168 
Davies, R. L., Sadler, E. M., & Peletier, R. F. 1993, 

MNRAS, 262, 650 
Davies, R. L. et al. 2001, ApJL, 548, L33 
Eisenstcin, D. J. ct al. 2003, ApJ, 585, 694 
Elson, R. A. W. 1997, MNRAS, 286, 771 
Faber, S. M. 1972, A&A, 20, 361 

Forbes, D. A., Brodie, J. P., & Larsen, S. ;. S. 2001, ApJL, 
556, L83 

Franx, M. & Illingworth, G. D. 1988, ApJL, 327, L55 
de Freitas Pacheco, J. A. & Barbuy, B. 1995, A&A, 302, 
718 

Fulbright, J. P. 2002, AJ, 123, 404 

Gonzalez, J. J. 1993, Ph.D. Thesis, University of 

California, Santa Cruz 
Goudfrooij, P., Hansen, L., Jorgensen, H. E., & Norgaard- 

Nielsen, H. U. 1994, A&AS, 105, 341 
Goudfrooij, P. & Emsellem, E. 1996, A&A, 306, L45 
Goudfrooij, P., Alonso, M. V., Maraston, C, & Minniti, 

D. 2001, MNRAS, 328, 237 
Goudfrooij, P., Strader, J., Brenneman, L., Kisslcr-Patig, 

M., Minniti, D., & Edwin Huizinga, J. 2003, MNRAS, 

343, 665 



16 



Puzia et al.: VLT Spectroscopy of Globular Cluster Systems 



Gregg,M. D., Ferguson,!!. C, Minniti,D., Tanvir,N.,& 

Catchpole,R. 2004, AJ, 127, 1441 
Greggio, L. 1997, MNRAS, 285, 151 
Greggio, L. & Renzini, A. 1990, ApJ, 364, 35 
Grillmair, C. J., et al. 1996, AJ, 112, 1975 
H arris, W .E. 1996 , AJ, 112, 1487, for the 1999 update see 

http : / /physun. physics .mcmaster . ca/ ~ harris / mwgc . dat 
Harris, W. E. 2001, Star Clusters, 223 
Harris, W. E. & van den Bergh, S. 1981, AJ, 86, 1627 
Harris, W. E. & Harris, G. L. H. 2002, AJ, 123, 3108 
Homeier, N., Gallagher, J. S., & Pasquali, A. 2002, A&A, 

391, 857 

Ibata, R. A., Gilmore, G., & Irwin, M. J. 1995, MNRAS, 
277, 781 

Jedrzejewski, R. & Schechter, P. L. 1988, ApJL, 330, L87 

Johnson, K. E., Vacca, W. D., Leitherer, C, Conti, P. S., 
& Lipscy, S. J. 1999, AJ, 117, 1708 

Kauffmann, G., White, S. D. M., & Guiderdoni, B. 1993, 
MNRAS, 264, 201 

Kauffmann, G. & Haehnelt, M. 2000, MNRAS, 311, 576 

Kissler-Patig, M., Richtler, T., Storm, J., & della Valle, 
M. 1997, A&A, 327, 503 

Kissler-Patig, M., Brodie, J. P., Schroder, L. L., et al. 
1998, AJ, 115, 105 

Khochfar, S. & Burkert, A. 2003, ApJL, 597, 117 

Koo, D. C. et al. 1996, ApJ, 469, 535 

Korn, A., Maraston, C, & Thomas, D. 2005, A&A sub- 
mitted 

Kraft, R. P. 1994, PASP, 106, 553 

Kuntschner, H. 2000, MNRAS, 315, 184 

Kuntschner, H., Smith, R. J., Colless, M., Davies, R. L., 

Kaldare, R., & Vazdekis, A. 2002a, MNRAS, 337, 172 
Kuntschner, H., Ziegler, B. L., Sharpies, R. M., Worthey, 

G., & Fricke, K. J. 2002b, A&A, 395, 761 
Larsen, S. S. & Richtler, T. 2000, A&A, 354, 836 
Larsen, S. S., Brodie, J. P., Beasley, M. A., & Forbes, 

D. A. 2002, AJ, 124, 828 
Larsen, S. S., Brodie, J. P., Beasley, M. A., Forbes, D.A., 

Kissler-Patig M., et al. 2003, ApJ, 585, 767 
Larson, R. B. 1975, MNRAS, 173, 671 
Lee, Y., Demarque, P., & Zinn, R. 1994, ApJ, 423, 248 
Lee, H., Yoon, S., & Lee, Y. 2000, AJ, 120, 998 
Macchetto, F., Pastoriza, M., Caon, N., et al. 1996, A&AS, 

120, 463 

Maraston, C. 1998, MNRAS, 300, 872 

Maraston, C. 2003, Extragalactic Globular Cluster 

Systems, Proceedings of the ESO Workshop held in 

Garching, Germany, 27-30 August 2002, p. 237 

Maraston, C. 2005, MNRAS accepted, ( |astro-ph/ 0410207) 
Maraston, C. & Thomas, D. 2000, ApJ, 541, 126 
Maraston, C, Kissler-Patig, M., Brodie, J. P., Barmby, P., 

& Huchra, J. P. 2001, A&A, 370, 176 
Maraston, C, Greggio, L., Renzini, A., Ortolani, S., 

Saglia, R. P., Puzia, T. H., & Kissler-Patig, M. 2003, 

A&A, 400, 823 
Maraston, C, Bastian, N., Saglia, R.P., Kissler-Patig, M., 

Schweizer, F. & Goudfrooij, P. 2004, A&A, 416, 467 
Matteucci, F. 1994, A&A, 288, 57 



McLaughlin, D. E. 1999, AJ, 117, 2398 

McWilliam, A. 1997, ARA&A, 35, 503 

O'Conncll, R. W. 1976, ApJ, 206, 370 

Osterbrook, D. E. 1989, Astrophysics of Gaseous Nebulae 

and Active Galactic Nuclei, University Science Books, 

Mill Valley 

Peebles, P. J. E. 2002, in A New Era in Cosmology, eds. 
Metcalfe, N. & Shanks, T., ASP Conference Series, in 
press 

Pettini, M. 2003, in Lectures of the XIII Canary Islands 

Winter School of Astrophysics, Cambridge University 

Press, ( |astro -ph/0303272) 
Phillips, M. M., Jenkins, C. R., Dopita, M. A., Sadler, 

E. M., & Binettc, L. 1986, AJ, 91, 1062 
Poggianti, B. M. et al. 2001, ApJ, 563, 118 
Puzia, T. H., Kissler-Patig, M., Brodie, J. P., & Huchra, 

J. P. 2000, AJ 120, 777 
Puzia, T. H., Saglia, R. P., Kissler-Patig, M., Maraston, 

C, Greggio, L., Renzini, A., & Ortolani, S. 2002, A&A, 

395, 45 

Puzia, T. H., Kissler-Patig, M., Thomas, D., Maraston, 
C, Saglia, R. P., Bender, R., Richtler, T., Goudfrooij, 
P., & Hempel, M. 2004, A&A, 415, 123, (Paper I) 

Puzia, T. H., Perrett, K. M., & Bridges, T. J. 2005, A&A, 
434, 909 

Reader J. & Corliss Ch.H. 1981, CRC Handbook of 

Chemistry and Physics 
Rich, R. M. et al. 1997, ApJL, 484, L25 
Rich, R. M. et al. 2003, AJ, in press 
Rosen, B. 1952, Constantes selectionnees; atlas des 

longueurs d'onde caracteristiques des bandes d'emission 

et d'absorption des molecules diatomiques, Hermann & 

C IC , Depositaries, Paris 1952 
Saglia, R. P., Maraston, C, Thomas, D., Bender, R., & 

Colless, M. 2002, ApJL, 579, L13 
Schweizer, F. 1997, ASP Conf. Ser. 116: The Nature of 

Elliptical Galaxies; 2nd Stromlo Symposium, 447 
Silk, J. 1977, ApJ, 214, 152 

Somerville, R. S., Primack, J. R., & Faber, S. M. 2001, 

MNRAS, 320, 504 
Surma, P. & Bender, R. 1995, A&A, 298, 405 
Sweigart, A. V. & Catelan, M. 1998, ApJL, 501, L63 
Thomas, D., Greggio, L., & Bender, R. 1999, MNRAS, 

302, 537 

Thomas, D., Maraston, C, & Bender, R. 2003, MNRAS, 
339, 897 

Thomas, D., Maraston, C, & Korn, A. 2004, MNRAS, 
351, L19 

Thomas, D., Maraston, C, Bender, R., & de Oliveira, 

C. M. 2005, ApJ, 621, 673 
Tinsley, B. M. 1972, ApJ, 178, 319 
Tinsley, B. M. 1979, ApJ, 229, 1046 

Tonry, J. L., Dressier, A., Blakeslee, J. P., et al. 2001, ApJ, 
546, 681 

Trager, S. C, Worthey, G., Faber, S. M., Burstein, D., & 

Gonzalez, J. J. 1998, ApJS, 116, 1 
Trager, S. C, Faber, S. M., Worthey, G., & Gonzalez, J. J. 

2000a, AJ, 119, 1645 



Puzia et al.: VLT Spectroscopy of Globular Cluster Systems 



17 



Trager, S. C, Faber, S. M., Worthey, G., & Gonzalez, 

J. J. ;. 2000b, AJ, 120, 165 
Trcu, T., Stiavelli, M., Casertano, S., M0ller, P., & Bcrtin, 

G. 1999, MNRAS, 308, 1037 
Tripicco, M. J. & Bell, R. A. 1995, AJ, 110, 3035 
Tully, R. B. 1988, Cambridge and New York, Cambridge 

University Press 
de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G. 

Jr., Buta, R., Paturel, G., and Fouque, P. 1991, Third 

Reference Catalogue of bright galaxies, Berlin: Springer 
White, S. D. M. & Rees, M. J. 1978, MNRAS, 183, 341 
Whitmorc, B. C. & Schweizer, F. 1995, AJ, 109, 960 
Worthey, G. 1994, ApJS, 95, 107 

Worthey, G., Faber, S. M., & Gonzalez, J. J. 1992, ApJ, 
398, 69 

Worthey, G., Faber, S. M., Gonzalez, J. J., & Burstein, D. 

1994, ApJS, 94, 687 
Worthey, G., & Ottaviani D.L. 1997, ApJS, 111, 377 



